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THE KINEMATIC PROPERTIES OF THE EXTENDED DISKS OF SPIRAL GALAXIES: A SAMPLE OF 

EDGE-ON GALAXIES 

Daniel Christlein 1-2 & Dennis Zaritsky 3 

Accepted for publication in the Astrophysical Journal 

ABSTRACT 

We present a kinematic study of the outer regions (R25 < R < 2 x R25) of 17 edge-on disk galaxies. 
Using deep long-slit spectroscopy (flux sensitivity ~ few 10~ 19 erg s _1 cm -2 arcsec -2 ), we search 
for Ha emission, which must be emitted at these flux levels by any accumulation of hydrogen due to 
the presence of the extragalactic UV background and any other, local source of UV flux. We present 
results from the individual galaxy spectra and a stacked composite. We detect Ha in many cases well 
beyond R25 and sometimes as far as 2 x i?25- The combination of sensitivity, spatial rcsultion, and 
kinematic resolution of this technique thus provides a powerful complement to 21-cm observations. 
Kinematics in the outer disk are generally disk-like (flat rotation curves, small velocity dispersions) at 
all radii, and there is no evidence for a change in the velocity dispersion with radius. We place strong 
limits, few percent, on the existence of counter-rotating gas out to 1.5i?25- These results suggest 
that thin disks extend well beyond R25', however, we also find a few puzzling anomalies. In ESO 
323-G033 we find two emission regions that have velocities close to the systemic velocity rather than 
the expected rotation velocity. These low relative velocities are unlikely to be simply due to projection 
effects and so suggest that these regions are not on disk-plane, circular orbits. In MCG-01-31-002 we 
find emission from gas with a large velocity dispersion that is co-rotating with the inner disk. 

Subject headings: galaxies: kinematics and dynamics - galaxies: formation - galaxies: evolution - 
galaxies: spiral - galaxies: structure 



1. INTRODUCTION 

Recent findings promulgate the notion that the 
baryonic component of spiral galaxies is significantly 
larger than suggested by classical size indicators such 
as the -B-band R25 radius (the isophotal radius where 
the surface brightness drops to 25 mag arcsec -2 ). In 
an era where the study of dark matter halos has shown 
spiral galaxies to be at least ten tim e s larger than 



the o ptical disk ( Zaritsky et al.l 119931 : iMcKav et all 
l2002t iPrada et al.l I2003D . one may not be surprised 



that baryons are also found at large radii, and radio 
observations of the 21-cm line of neutral hydro- 
gen have anticipated this suspicion by many decade s 
(|van de Hulst. Raimond. k Woerde^ll957t[Dieterill962l ). 
More recently, striking new discoveries are showing that 
these baryons, which are found in all three principal 
"phase s" (stars (iFerguson et al.ll2002l 120051: iPohlen et all 
2002t [Zibetti et all 12004 IZibetti fc Ferguson! 12004 



Bland-Hawthorn et all l2005t lEUis k Bland-Hawthornl 
20061). gas dPhookun. Voeel k Mundvl 119931 : 
Bland -Hawthorn et all Il997t iFerguson et all Il998allbl: 
Thilker et al.ll2005h. and dust (jNelson. Zaritsky k Cutril 
1998t lAlton et alJ 1200 ll )) are organized in such a way 
that they are not merely the flotsam left as the forming 
galaxy falls into the dark matter potential well, but 
rather the outer extent of the structures manifested at 
smaller radii. How these structures transition from the 
baryonic disk, where the baryonic mass has a sharply 
peaked distribution of specific angular momenta, to the 
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halo of the galaxy, where the measurable baryons, those 
in globular clusters and satellite galaxies, have a broad 
distribution of specific angular momenta, is still unclear. 

These two realms are not isolated. Hierarchical growth 
models predict that the galaxy grows by accreting new 
material from the halo and beyond. As the infalling outer 
material crashes onto the outer banks of a disk galaxy, 
it should leave kinematic signatures in the form of lower 
rotational velocities, increased velocity dispersion, and 
possibly even counte rrotation in the inter f ace re gion be- 
tween disk and halo. Ivan den Bosch et al.l (|2004l ) predict 
that, upon the formation of a galaxy, up to 40% of the 
gas in a halo has negative specific angular momentum. If 
such gas stands any chance of survival, it would be in the 
outermost regions of late-type galaxies, where dynamical 
timescales are much larger than in the inner disk, where 
such signatures in the gas kinematics would be quickly 
wiped out. 

The outer disk may also be a significant repository of 
positive angular momentum. If the surface mass density 
of the disk follows an exponential law and the rotation 
curve remains flat, then around the radius -R25, the to- 
tal angular momentum will scale roughly in proportion 
to the radius. This additional angular momentum in 
the outer disk sets an even higher bar for hierarchical 
growth models, which underpredict the an gular momen- 
tum of galaxies (jvan den Bosch et al.ll2004T ). A kinematic 
inventory of the outskirts of disks is therefore an impor- 
tant aspect of a full accounting of the angular momentum 
of disk galaxies. 

Although radio observations have been the stalwart of 
this field for over three decades, and have confirmed that 
many galaxies are embedded in extended disks of neutral 
hydrogen, they have some potentially important limita- 
The large beam size of radio telescopes means 
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that, at each pointing, material from a wide range of 
radii, elongations, and velocities contributes to the ob- 
served signal, making detailed modeling of the emission 
line necessary. If gas in the outskirts of galaxies is dis- 
tributed in clumps rather than smoothly, such detailed 
modeling becomes very difficult, and small sources with 
anomalous kinematics may be indistinguishable from the 
velocity distribution that would be expected across the 
beam anyway. In addition, if t he gas is ionized due to ir- 
radiation by the inner galaxy (|Bland-Hawthornlll998f) or 
the intergalactic ultraviolet background (|Sunvaevlll969t 
lMalonevlll99l : iCorbelli fc Salucc1ll993h then there will 
be no 21-cm emission. Furthermore, radio interfero- 
metric observations require substantial observational ef- 
forts to achieve favorable combinations of velocity reso- 
lution, spatial resolution, and high sensitivity, and usu- 
ally strike a trade-off particularly b etween spatial reso- 
lution and sensitivity (for example, iMartinl (|1998l ) lists 
no arcsecond-scale 21-cm maps of extragalactic objects). 
Even so, interferometric 21-cm observations for even very 
nearby galaxies often extend over hours or even days. Fi- 
nally, 21-cm line observations provide no ancillary infor- 
mation, i.e., other emission lines, that could be used for 
line diagnostics of density and metallicity. 

An underexplored way to study the kinematics of the 
outer regions of gaseous disks is by detecting optical 
emission, for example, the Ha line. Due to the low 
level of star formation i n the outer disks of galaxies, typ- 
ical H a rotation curves (jVogt et al.ll2004t ICatinella et alJ 
2006) end well within the optical radius R25 of the 
disk. Nonetheless, Ha emission is expected at some level 
from any accumulation of hydrogen due to the ubiq- 
uity of ionizing flux throughout the universe not only 
from local star formation, but also from nearby galaxies, 
AGN, and the metagalactic UV background, even though 
such emission may only occur at extremely faint sur- 
face brightness levels (~ 10 -19 erg s _1 cm -2 arcsec~ 2 ). 
First attempts to detect suc h emission were made by 
iBland-Hawthorn et al.1 (|1997l ) with a Fabry-Perot inter- 
ferometer, revealing outer-disk emission in NGC 253 even 
beyond the truncation radius of the neutral hydrogen 
disk at a level of a few 10~ 19 erg s _1 cm -2 arcsec~ 2 . 
With modern 6-8 m class telescopes and efficient detec- 
tors, it is possible to push the detection threshold even 
with long-slit spectrographs to sufficiently low sensitivi- 
ties to detect Ha well beyond the bright stellar disk, in a 
realm usually only probed by 21-cm observations, and at 
the same time make use of the more robust sky subtrac- 
tion that long-slit spectrography offers over Fabry-Perot 
observations. Long exposure times (several hours) are 
required for a spectroscopic detection of the outer disk, 
but if successful, such detections will provide high spa- 
tial and reasonably high kinematic resolution, and a host 
of ancillary information, such as a measure of the stellar 
continuum and other emission lines, most notably [SII] 
and [Nil]. 

Most searches for Ha emission from the outer disks in 
the past, including those listed above, have made use of 
narrow-band imaging. However, for line emission, opti- 
cal spectroscopy is capable of achieving a much higher 
S/N per pixel than imaging, as it allows us to isolate 
the line flux from the sky and stellar continuum back- 
ground much more reliably. We can therefore probe 
much fainter structures in the outer disk, including not 



only selected bright HII regions, but also the diffuse in- 
terstellar medium. 

We report here on the first results of an observational 
campaign to study the outer disks of galaxies with long- 
slit spectroscopy. Specifically, in this paper, we discuss 
constraints on the kinematics of the outer disk. In §2.11 
we discuss what we might detect in Ha in the outer disks. 
In Sj2l we describe our observational campaign. In Sj3j we 
outline our reduction and analysis pipeline. Our results 
and discussions thereof are presented in [21 We examine 
rotation curves out to radii of 2 x R25, i.e., our best cases 
are competitive with the extent reached by radio inter- 
ferometric studies. We specifically search for kinematic 
anomalies, i.e., gas that is systematically deviating from 
co-rotation, taking advantage of the high spatial resolu- 
tion of optical spectroscopy. We also discuss what con- 
straints can be placed on the profile of a potential dark 
matter halo at large radii from these data. In addition, 
we use a stacked composite spectrum to examine the ve- 
locity dispersion of gas in the outer disk and to search 
for gas on counterrotating orbits. 

Subsequent papers will discuss the physical state of 
the outer-disk gas, including metallicity and density and 
constraints on the intergalactic UV background, and will 
examine the relation between Ha emission and stellar 
populations in the outer disk. 

2. OBSERVATIONS 

2.1. Sources of Ha emission in the outer disk 

The brightest and most easily detected source of Ha 
emission are local star formation regions in the outer 
disk. Such regi ons have been i dentifie d using narrow- 
band imaging (jFerguson et al.l 1998al) and examined 
spectroscopically (jFerguson et alllT998bl ). The luminosi- 
ties of these objects are in the range of 0.1 — 8 x 10 38 
erg s _1 , although the lower limit is observational rather 
than physical. Using a long-slit spectrograph on an 8-m 
class telescope for a galaxy at a cosmological redshift of 
cz=4,000 km s _1 , such star formation regions are typi- 
cally detectable within several minutes. However, these 
objects highlight only small, selected regions, and may 
not cover the full extent of the gas distribution. Further- 
more, in a blind survey with a long-slit spectrograph, it 
is a matter of chance whether they fall onto the slit. 

A more comprehensive study of the baryon budget 
is possible if the outer disk is irradiated uniformly by 
a sufficiently strong exterior ionizing flux. The inner, 
star-forming disk itself can constitute such a source, es- 
pecia lly if the disk is slightly warped (Bland-Hawthorn 
1998). The Ha surface brightness in this case depends 
strongly on the star formation rate in the inner disk, the 
UV escape fraction, and the geometry of the disk, so 
the necessary exposure times can vary between minutes 
and several hours. A second such source is the inter- 
galactic UV background. The intensity of this radiation 
field is still uncertain; a theoretical predictio n estimates 



1.3+" 8 x 10- 23 



iw erg s x cm " sr x Hz 1 (|Shull et al.l 
1993). IVogel et all (|1995f ) provide an upper limit of 

1 Hz -1 , and a recent mea- 



< 8 x 10 23 erg s 1 cm 2 sr 



surement from (Scott et alJ (|2002f ) finds 7.6l^ x 10~ 23 



erg s 1 cm 2 sr 



1 Hz -1 . Ionization by the intergalactic 
field sets the lowest surface brightness limit of Ha emis- 
sion, but even it is within the reach of large telescopes 
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and modern instruments. Because Ha emission stimu- 
lated by the UV background provides the most uniform 
and unbiased probe of baryonic material — wherever 
there is hydrogen, there must be some Ha emission — , we 
set our exposure times to provide a realistic chance of de- 
tecting background-ionized hydrogen for at least some of 
our target galaxies. The Ha surface b rightness level cor - 
responding to the background flux of lScott et al.1 (|2002D 
is ~ 2^_\ x 10 -19 erg s _1 cm -2 arcsec -2 (for deriva- 
tions, see, e.g., iMalonevI (|1993h : iBland-Hawthorn et al.1 
(1997)). This sensitivity is attainable with an 8-m class 
telescope within several hours of exposure time. 

2.2. Target Selection 

We select target galaxies using the following criteria: 

1) visibility — targets should be observable during the 
entire night and have the lowest possible mean airmass, 

2) redshift — targets should have recessional velocities 
between 3,600 and 12,000 km s _1 for which the redshiftcd 
Ha line falls into a dark and flat region of the night sky 
spectrum, 3) angular size — the slit length should cover 
the region of interest, ideally the entire galaxy, and leave 
room for blank sky, 4) morphology — targets should be 
predominantly late-type galaxies (Sc-Sd) because these 
are generally gas-rich, 5) H I content — targets are pre- 
ferred if they have existing detections in the literature, 
6) inclination — targets must be close to edge-on con- 
figuration to maximize the chance that disk regions will 
lie within the slit, and 7) environment — targets must 
be relatively isolated both in regards to possible galaxy 
companions, but also in terms of bright stars that would 
affect the observations. 

These criteria are somewhat flexible, and in particular 
we do include an Sb, an Sab, and a gas-rich SO (which 
will be discussed in a separate paper; Kannappan et al., 
in prep.). Our sample does not constitute a complete or 
systematic sample in any way, but we emphasize that, 
with the exception of the SO, all of them were selected 
for appearing as undisturbed as possible. In practice, the 
axis ratios (b/a) of our target galaxies are < 0.2, with 
the exception of one object — UGC 10972, for which it 
is 0.27. The inclination angle in Tabled] is calculated di- 
rectly from the published axis ratio, neglecting the finite 
thickness of the disk, and is typically around 80 degrees. 
If more realistic values for the thickness are adopted (i.e., 
c/a < 0.2) all galaxies with the exception of UGC 10972 
are consistent with exact edge-on configuration. This 
becomes an important distinction, because in a perfectly 
edge-on configuration, the line of sight traverses the en- 
tire disk and thus a wide range of line-of-sight velocities, 
leading to possible ambiguities in the interpretation of 
the kinematic measurement. In a slightly inclined sys- 
tem, the slit probes a much smaller range of depths. We 
will discuss this topic more extensively below in connec- 
tion with measurement uncertainties and biases. 

The selected galaxies are listed in Table [TJ along with 
their redshift, the isophotal radius R25, the inclination 
angle, morphological type, and approximate absolute B- 
band magnitude, as gleaned from the NASA-IPAC Ex- 
tragalactic Database (NED) and calculated assuming 
a ACDM cosmology with Ho=71 km s _1 Mpc -1 and 
fi ro =0.27, S1 A =0.73, and Q k =0. 

2.3. Instruments 



Our observational campaign uses long slit spectro- 
graphs. Their advantages for our program include an 
extended field of view that covers the entire major axis 
of a galaxy simultaneously plus a substantial sky back- 
ground, and sufficient velocity resolution to determine 
rotation curves. 

A moderate dispersion spectrograph (of the order 
lA/pixel) is suitable for this project. High dispersion 
spectroscopy will improve the kinematic precision, but 
dilute the Ha line and so increase the relative contribu- 
tion of read noise. Low resolution spectroscopy would 
decrease the kinematic precision below what is necessary 
to measure kinematics internal to the galaxy; further- 
more, since the width of the line image on the detector 
is set primarily by the slit width, low-resolution spec- 
troscopy would increase the contribution of sky across 
this width. Anamorphic demagnification, which exists 
in some of the spectrographs used for our program (the 
B&C and GMOS-S spectrographs) further aids our ob- 
servations by reducing the effective slit image width, and 
thus minimizes noise contributions from read noise and 
sky background. 

We observed using the following instruments: a) the 
Boiler & Chivens spectrograph on the Magellan I (Baade) 
telescope using the 1200 1/mm grating blazed at 5500 A 
with a 1 arcsec slit on the nights of 26-28 Mar 2003 and 
16-17 Apr 2004, b) the Red Channel Spectrograph on the 
MMT using the 1200 1/mm grating blazed at 7700 A with 
a 1 arcsec slit on the nights of 21-22 May 2004, c) FORS-1 
on the VLT-2 with a 1 arcsec slit and the 600R grism with 
a dispersion of 1.038 Apixer 1 (unbinned) on Nov 16-18, 
2004, and d) GMOS on Gemini-South with a 1 arcsec 
slit and the R831 grating with a dispersion of 0.34 A per 
unbinned pixel on Apr 6-8, 2005. Aside from the throug- 
put and sensitivity differences, the principal difference in 
the setup between the GMOS and FORS1 observations 
on the one hand and the B&C and Red Channel obser- 
vations on the other hand lies in the smaller slit length 
available on the latter two spectrographs. In these cases, 
we first obtained a short (~ 5 — 10 min) spectrum with 
the galaxy centered on the slit, followed by much longer 
exposures with an offset along the major axis, usually in 
only one direction. 

Individual exposures, typically 20-30 min in length on 
the MMT and Magellan, and one hour on the VLT and 
Gemini, are combined to result in the total observing 
time given in Tabled Additional details of the observing 
runs are also provided in Table [5] 

3. DATA ANALYSIS 

3.1. Reduction 

We reduce the long-slit spectra from the MMT, Mag- 
ellan, and the VLT with the standard IRAF reduction 
tasks. For Gemini data, we use the GEMINI package in 
version 1.7 4 For the Gemini and VLT data, we determine 
wavelength solutions from arc spectra obtained immedi- 
ately before or after the science exposures. If there are 
any indications that the wavelength solution changed be- 

4 During the reduction, we discovered an error in the GSFLAT 
algorithm that caused the task to carry out the flatfielding with 
an individual flatfield, instead of a combined master flatfield; this 
problem was quickly patched by the Gemini observatory staff, and 
we continued the reduction with the patched version of the package. 
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Table 1. Basic Galaxy Data 



Galaxy 


cz 


K-25 


K-25 


b/a 


i 


Type 


Mb 




n — 1 1 
[km s 1 ] 


[arcsec] 


[kpc] 




[°] 






ESO 201-0022 


4068 


75.7 


20.7 


0.14 


82 


Shr HIT 

JUL J. ±±± 


— 19.0 


ESO 299-O018 


4821 


61.5 


19.9 


0.11 


84 


Sc 


— 18.4 


ESO 478 ooi 1 


5216 


42.1 


14.7 


0.12 


83 


Scd 


— 18.1 


TP 20^8 


1379 


93.1 


8.6 


0.15 


81 


Sc 


— 17.5 




2348 


58.7 


9.3 


0.16 


81 


Sd 


— 17.7 


eso 380-no23 


2753 


41.6 


7.7 


0.18 


80 


Scd 


— 17.3 


f,so 385-ooos 

UuW JUL) vJUUO 


3756 


62.2 


15.7 


0.10 


84 


Sd 


— 18.3 


IC 4393 


2753 


80.2 


14.9 


0.11 


84 


Scd 


-18.5 


MCG -01-31-002 


5741 


42.4 


16.2 


0.17 


80 


Scd 


-19.6 


ESO 445-G081 


4352 


60.8 


17.8 


0.18 


80 


Sbc 


-19.8 


ESO 445-G085 


4285 


57.4 


17.0 


0.06 


86 


Scd 


-17.7 


ESO 568-G010 


5507 


52.4 


19.3 


0.06 


87 


Sd 


-18.1 


ESO 381-G045 


7129 


33.0 


15.6 


0.20 


78 


Sb 


-19.6 


UGC 09138 


4601 


58.5 


18.1 


0.12 


83 


Sc 


-18.9 


UGC 09780 


5173 


62.7 


21.7 


0.12 


83 


Scd 


-19.0 


UGC 10453 


4352 


41.4 


12.1 


0.10 


84 


Sd 


-18.2 


UGC 10972 


4652 


75.3 


23.5 


0.27 


74 


Scd 


-20.0 



No te. — Data are drawn from NED. Isophotal radii R25 in the _B-band are taken from the ESO-LV catalog (|Lauberts fc Valentiinl 
[19891) for all ESO and IC galaxies, and from the RC3 D25 catalog ide Vaucouleurs et al.|[l99ll) for MCG -01-31-002 and the UGC galaxies. 
Inclination angles are calculated assuming the idealized case of an infinitesimally thin disk. 



Table 2. Summary of Observing Runs 



Galaxy 


Offset 


Telescope 


Date 


iEXP 
[sec] 


ESO 201-G022 


central 


VLT/FORS1 


Nov 2004 


37,329 


IC 2058 


central 


VLT/FORS1 


Nov 2004 


12,399 


ESO 299-G018 


central 


VLT/FORS1 


Nov 2004 


12,199 


ESO 478-G011 


central 


VLT/FORS1 


Nov 2004 


10,300 


ESO 323-G033 


central 


Gemini-S/GMOS 


Apr 2005 


39,600 


IC 4393 


central 


Gemini-S/GMOS 


Apr 2005 


15,400 


ESO 380-G023 


central 


Gemini-S/GMOS 


Apr 2005 


10,800 


ESO 385-G008 


central 


Gemini-S/GMOS 


Apr 2005 


7,200 


ESO 381-G045 


NE 


MAG/B&C 


Mar 2003 


50,400 


ESO 385-G008 


NE 


MAG/B&C 


Apr 2004 


4,800 


ESO 445-G081 


NE 


MAG/B&C 


Apr 2004 


19,200 


ESO 445-G085 


NE 


MAG/B&C 


Apr 2004 


18,600 


ESO 568-G010 


NE 


MAG/B&C 


Apr 2004 


15,800 


IC 4469 


NW 


MMT/Red 


May 2004 


8,284 


MCG-01-31-002 


NE 


MAG/B&C 


Apr 2004 


9,000 


UGC 09138 


NW 


MMT/Red 


May 2005 


3,600 


UGC 09780 


SE 


MMT/Red 


May 2004 


3,600 


UGC 09780 


NW 


MMT/Red 


May 2004 


4,800 


UGC 10453 


SW 


MMT/Red 


May 2005 


4,800 


UGC 10972 


NE 


MMT/Red 


May 2005 


4,800 


UGC 10972 


SW 


MMT/Red 


May 2005 


3,600 



Note. — "Offset" indicates where the slit center was positioned 
relative to the galaxy (offsets are always along the major axis). 



tween exposures, we use separate solutions for different 
images or groups of images. For the Magellan and MMT 
spectra, we use sky lines to wavelength calibrate. We 
sky subtract using the sky spectrum from the ends of 
the slit. The region that we consider to be devoid of 
emission from the galaxy itself is determined by eye, and 
typically corresponds to a fraction of at least half the 
slit length. We continuum subtract by fitting low-order 
polynomial functions along the dispersion axis. 

3.2. Measurements 



We use the Ha line for the measurements of the gas 
kinematics in the outer disk. Although other lines ([Nil], 
[SII]) are observed, these are usually very weak in the 
outer regions of our target galaxies. To detect and ex- 
tract the Ha line, we begin by applying a smoothing filter 
with a width of ~ 1.5 arcsec in the spatial direction to 
the spectrum. This scale is about twice the typical see- 
ing and is chosen as a compromise to detect both small- 
scale features and attain sufficient sensitivity. Starting 
at a local maximum near the wavelength range where Ha 
emission is expected, we define a search window whose 
width along the dispersion axis narrowly encloses the 
entire line, and determine the mean intensity-weighted 
wavelength within this window as the new estimate for 
the line centroid. We iterate this process, setting the new 
estimate for the centroid as the new center of the search 
window. This process is carried out interactively, us- 
ing software that provides a direct visual feedback of the 
line extraction process and allows us to correct errors. 
Failures in the automatic centroid search are corrected 
manually by restarting the iteration near the correct po- 
sition. At very large radii, where visual inspection often 
yields no significant peak, the decision where to start the 
iteration admittedly becomes arbitrary, but in such cases 
the resulting detections generally turn out to be below 
the significance threshold as well; i.e., the risk of selecting 
spurious peaks due to "hunting for significance" , is not a 
major problem in practice. After convergence, fluxes are 
determined by integrating the emission within the search 
window. 

A source of systematic error in the determination of 
fluxes, and thus of the significance of a detection, is broad 
Ha absorption by the underlying stellar population. In 
the bright inner disks, Ha absorption is negligible, com- 
pared to the Ha emission, but in some cases towards 
the outer edge of the disk, where the stellar continuum 
is noticeably more extended than Ha emission, it can 
lead to a systematic underestimate of the Ha flux, or 
even suppress the Ha emission completely. To remove 
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this local background, we define regions along the disper- 
sion axis between the Ha emission line and the [Nil] 6548 
and [Nil] 6584 lines, where the Ha absorption is observed, 
at fixed wavelength offsets relative to the Ha centroid. 
We then fit — for each row of the input spectrum - 
a Gaussian to the Ha absorption trough in these local 
background regions. Typically, the Gaussian parameters 
(width, amplitude, centroid) exhibit large scatter from 
one row to the next. In a next step, we therefore fit 
the parameters of the individual Gaussian fits as a func- 
tion of the position along the slit with simple low-order 
functions (an arctan function for the centroid, Gaussian 
for the width, and Gaussian or exponential, as suggested 
by the data, for the amplitude). With the parameters 
interpolated from these fits, we model the local Ha ab- 
sorption trough in each row. We also review the quality 
of the local background model and apply manual correc- 
tions on a row-by-row basis where necessary. In particu- 
lar, we adopt a conservative approach towards outlying 
emission regions, and in general do not accept an iso- 
lated detection at large distances from the galaxy as real 
if it is pushed over the significance threshold only by the 
local background subtraction (unless the presence of Ha 
absorption is unmistakable). 

3.3. Uncertainties 

Uncertainties in the measured flux are determined in 
two ways: First, we apply a bootstrapping algorithm, 
which randomply resamples the CCD rows within the 
area covered by our measurement window, and measures 
the flux from this random sample. By repeating this 
and determining the standard deviation of the results, 
we obtain a first estimate for the flux uncertainty. While 
this algorithm accounts well for variations on pixel-to- 
pixel scales, including Poisson errors in the Ha flux, it is 
insensitive to fluctuations in the sky or continuum back- 
ground that occur on scales comparable to that of the 
measurement window or larger. Such fluctuations fre- 
quently result from flatfielding imperfections and ulti- 
mately a more severe limitation of our ability to detect 
faint Ha emission than the Poisson errors in the sky back- 
ground. In order to quantify the uncertainties in the 
sky background, we define several control regions along 
the dispersion axis that do not contain prominent emis- 
sion or absorption lines. We then slide our measurement 
window across this background region, determining the 
integrated flux at each position. Since these control re- 
gions are typically at least an order of magnitude larger 
than the measurement window, the standard deviation 
of the measured residual fluxes around their mean is a 
good measure of the typical errors in the sky background 
within the measurement window. As our final flux uncer- 
tainty, we adopt the larger of the two estimates described 
in this paragraph. Usually, the uncertainties determined 
after the first method are larger in the central parts of 
the galaxy, where the Poisson errors in the Ha emission 
dominate, while at large radii, the large-scale fluctua- 
tions in the sky background are the dominating source 
of error. 

In Fig. [TJ we show how the flux uncertainty varies as 
a function of radius, using one of our deepest and most 
extended spectra (ESO 20I-G022) as an example. Bold 
data points mark flux measurements that we consider 
significant, and small dots denote insignificant measure- 




Fig. 1. — Comparison of the two uncertainty estimates in the 
case of ESO 201-G022. The horizontal axis shows the distance 
from the galaxy center in units of i?25 • The vertical axis shows the 
Ha surface brightness and is linear below a value of 1, and decad- 
ically logarithmic above (i.e., an abscissa of 2 signifies a value of 
10 etc.). Bold data points show significant flux detections, small 
dots insignificant ones. The dotted line indicates flux uncertainties 
from bootstrapping, and the solid line shows background uncertain- 
ties. In the inner disk, uncertainties based on the bootstrapping 
method, which include Poisson errors from Ha emission and the 
continuum, dominate, but the outer disk is dominated by errors in 
the background, which includes not only Poisson noise in the sky 
background, but also correlated errors, such as flatfielding imper- 
fections. We use the larger of the two to evaluate the significance 
of a detection. 



ments. The dotted, thin line at the bottom represents 
the flux uncertainties as determined by the first method 
(bootstrapping within the measurement window), while 
the bolder, continuous line shows the background flux 
uncertainties, obtained by the second method described 
above. The plot is linear below an ordinate of 1, and 
decadically logarithmic above (i.e., a value of 1 repre- 
sents a measurement of I of the measured units, while a 
value of 2 represents 10 units, and a value of 0.5 repre- 
sents 0.5 units). 

For each spectrum, we determine a detection thresh- 
old by fitting a low-order functional form to the extracted 
data points and then plotting the residual velocity ver- 
sus the fit as a function of the significance of the mea- 
sured flux. In order to maintain a conservative approach, 
this detection threshold may vary from spectrum to spec- 
trum, depending on the amount of flatfielding imperfec- 
tions, skyline residuals, and other factors, and is between 
3 and 7 a for our sample. 

Uncertainties in the velocity centroid are determined 
analogously to the first method, i.e., we measure the 
standard deviation in the recovered line centroids as we 
apply a bootstrapping algorithm that randomly resam- 
ples the CCD rows within our measurement window (the 
intensity distribution along the dispersion axis is not af- 
fected by the bootstrapping) and estimates the disper- 
sion of the recovered centroids. 

This method reflects only the statistical uncertainty in 
the centroid position, not in the actual line-of-sight ve- 
locity. In practice, another source of error dominates: if 
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the slit is not uniformly or at least symmetrically illu- 
minated, the recovered line centroid will not reflect the 
actual line-of-sight velocity, and surface brightness fluc- 
tuations across the slit width may be misinterpreted as 
velocity shifts. With the spectrographs used in our study 
(dispersion and slit width), this effect can theoretically 
introduce deviations in the recovered velocity of several 
tens of km s _1 . Higher dispersion, narrow slits, and 
stronger anamorphic demagnification of the slit image, 
but also larger seeing, help reduce this effect. In general, 
however, this issue is only a problem in determining the 
velocities of individual emission regions. Due to the small 
angular width of the slit relative to the angular sizes of 
the galaxies, sources are distributed stochastically across 
the slit, so that, for large samples — i.e., those binning 
over a large spatial extent, or those co-adding multiple 
galaxies — , this effect will not introduce strong system- 
atic deviations in the rotation curve. We quantify the 
effect on the line centroid measurements of individual 
regions by examining the scatter of the measured line 
centroids about the fitted rotation curve. The measured 
scatter is on the order of 10-20 km s _1 . This result does 
not preclude the possibility that actual velocity fluctu- 
ations with such amplitudes may exist, but we do not 
attempt to distinguish between these two possibilities in 
this analysis; such a distinction would only be possible 
via the shape of the line profile, which is only viable in 
the high-surface brightness inner regions of the galaxy. 

In addition to this effect, which introduces a random 
uncertainty in the recovered line-of-sight velocity, there 
is a potential systematic source of error. Deriving the 
line-of-sight velocity from the intensity-weighted mean 
wavelength, as we have done, is an approximation. In 
reality, any line of sight traversing an extended part of a 
strongly inclinded disk will intercept parcels of emitting 
gas at a range of line-of-sight velocities, which are not dis- 
tributed symmetrically with respect to the rotational ve- 
locity at the projected radial distance. This usually leads 
to an underestimation of the true rotational velocity by 
generating a tail in the emission profile towards lower ve- 
locities. For geometric reasons, the effect is stronger near 
the center of the galaxy than on the outskirts (it depends 
primarily on the ratio of the slit width to the projected 
distance from the center of rotation). This effect is a 
serious concern in radio astronomy, where beam widths 
are large and velocity differentials across the beam can 
be significant. With optical spectroscopy, it is generally 
smaller by an order of magnitude, because the width of 
the spectrograph slit is small compared to the typical 
diameter of a galaxy. To quantify this bias, we have esti- 
mated the ratio of the recovered mean rotational velocity 
to the true rotational velocity, using simple models with 
a constant rotational velocity and gas surface density, as 
well as finite thickness of the gaseous disk. For an in- 
clination angle of 85°, the deviation between the mean 
flux- weighted velocity, which we recover, and the true 
rotational velocity is within our typical uncertainties as 
long as the ratio of the projected radial distance to the 
slit width is greater than 30. This is the case at the 
R25 radius for all galaxies in our sample. For all other 
galaxies, the effect becomes a concern at R25 only for 
inclination angles of 88° or more. This is a conservative 
estimate; in reality, the Ha flux drops steeply at large 
radii, so that the flux in the low-velocity tail is further 



suppressed. Final affirmation that this effect is no great 
concern (even though most of our sample is consistent 
with exact edge-on confirmation) comes from the data 
themselves: In our composite spectrum (which we dis- 
cuss in detail in §4.5. lj) . the Ha line appears symmetric 
which no indication of a strong low-velocity tail at any 
radius. 

In edge-on galaxies, dust is known to affect optical ro- 
tation curves: At small projected distances from the nu- 
cleus, dust extinction suppresses emission from the inner 
parts of the galaxy, increasing the relative contribution 
of emission from larger radii along the same line of sight, 
which tends to have a smaller line-of-sight velocity; in 
other words, the low- velocity tail discussed in the pre- 
vious paragraph is emphasized relative to the peak of 
emission near the rotational velocity. This is known to 
severely bias the slope of th e rotation curve at small radii 
(|Giovanelli &; Havnesll2002| ). For this reason, our survey 
is not ideally designed for determining the shape of the 
rotation curve at small radii. At larger radii, however, 
all these biasing effects are weaker: not only is the veloc- 
ity gradient within the field of view smaller, as argued 
above, but we also expect dust absorption to be a very 
minor problem around R25 and beyond, where the focus 
of our study lies (which is also noted by Giovanelli & 
Haynes) . 

An additional source of systematic errors in optical 
spectra arises from possible misalignment between the 
slit and the galaxy nucleus; however, this effect is only 
a concern at extremely small radii. Since the present in- 
vestigation does not focus on the shape of the rotation 
curve near the nucleus, we neglect this effect. 

Finally, if there are strong warps present in the outer 
disk and their line of nodes is aligned perpendicular to 
the line of sight, then the true inclination angle of the 
outer gaseous disk may be different from that derived 
from the optical axis ratio, and the observed line-of-sight 
velocity may be an underestimate of the true rotational 
velocity in the warped disk. This possibility cannot be 
ruled out from optical data along a one-dimensional slit 
alone. However, in our data, we do not see the signature 
of such an effect (a sudden drop in the recovered line- 
of-sight velocity), which may indicate that either such 
configurations (i.e., warp tilts on the order of 25° or more, 
required to introduce a significant dropoff in the velocity) 
are rare, or that our rotation curves do not extend to 
such radii. The possibility exists that an actual rise in 
the rotation curve is cancelled out exactly by this tilting 
effect, but it is highly unlikely that these would conspire 
to reproduce an undisturbed-looking rotation curve in 
every single case of our sample. 

3.3.1. Flux Calibration 

We flux-calibrate our spectra using the standard stars 
GD108 and LTT9491 for our VLT data from November 

2004, and LTT6248 for the Gemini-S data from April 

2005. We calculate fluxes by integrating over the entire 
line width, and multiplying with the calibration factor 
(which is measured in erg cm~ 2 count -1 ). We then con- 
vert the flux into a surface brightness, using the slit width 
(1 arcsec for all galaxies in this sample), equivalent expo- 
sure time, and platescale in the spatial direction. Flux 
calibrations from any one run are consistent to within 
5-10%. 
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For data from Magellan and the MMT, no spectropho- 
tometric standard star observations are available. We 
use the one object (ESO 385-G008) observed both with 
Gemini-South and Magellan, to derive an approximate 
relative calibration factor between the two optical sys- 
tems, and determine that factor to be ~ 4.11 (i.e., for the 
same incident flux, Gemini-S/GMOS measure a count 
rate 4.11 times higher). We assume the same value for 
the MMT (the diameters of the MMT and Magellan pri- 
mary mirrors are comparable; however, the Red Channel 
spectrograph on the MMT is likely to be more sensi- 
tive than the B&C spectrograph on Magellan, therefore 
this assumption is very approximate). We are aware of 
the inherent uncertainties in this relative calibration, but 
in this paper, we only use the absolute flux to assign a 
statistical weight to each spectrum when superimposing 
data from different runs and/or objects. Therefore, the 
actual impact of the calibration on our results is minimal. 

4. RESULTS 

4.1. Rotation Curves 

The rotation curve at large radii is a kinematic probe of 
the mass profile in a region where luminous matter con- 
tributes little to the total mass budget. Steep declines 
in the rotation curve would therefore indicate a break or 
truncation in the dark matter halo density profile. Al- 
though the general existence of massive dark matter ha- 
los at radii much larger than tho se probed here is inf erred 
by studies of satellite galaxie s (|Zaritskv et al.lll993t ) and 
weak lensing ( jBraineraf 19961 ) . these results are statistical 
in nature. Ha rotation curves, along with 21-cm rotation 
curves, are the sole avenue to study the halo mass profiles 
of a large sample of individual spiral galaxies, albeit at 
radii smaller than those probed by satellite galaxies and 
lensing. Furthermore, kinematic anomalies, such as an 
increased velocity dispersion, or individual emission re- 
gions deviating from co-rotation, are possible signatures 
of the hierarchical build-up of a galaxy predicted under 
the current cosmological paradigm. At large radii, where 
dynamical timescales are long, such signatures are pre- 
served for longer periods than in the inner disk. 

In this section, we present our results for the rotation 
curves of the individual galaxies in our sample, based 
on the measured Ha velocity centroids. Our sample 
is designed to constrain the rotation curve at relatively 
large radii, close to and beyo nd the isophotal rad ius R25. 
Stellar continuum emission (|Pohlen et al.ll20"02T) and, in 
fact , resolved stellar po pulations in local group galax- 
ies (jFerguson et al.ll2~002D have been measured to signifi- 
anctly fainter surface brightnesses, and R25 therefore cer- 
tainly does not represent a strict cut-off radius of the stel- 
lar disk. Nonetheless, this radius is easily reproducible 
and contains ~ 90% of the total flux of a galaxy. Typi- 
cally, i?25 is approximately twice as large as the half-light 
radius, and 5-10 times the surface brightness scale length 
in the i?-band. It is therefore a very suitable measure of 
the radius of the bright, stellar disk of a galaxy, and we 
cast our measurements of radial distances in units of R25 
to facilitate the comparison between different galaxies. 

The values of R25 used here are drawn from the 
NASA Extragalactic Database (NED) and come from 
two sources: the Surface Photometry Catalogue o f 
the ESO-Uppsala Galaxies (jLauberts fc Valentiinlll989l ) 
and the Third Reference Catalog of Bright Galaxies 



(|de Vaucouleurs et al.1 Il99l[ ). Table Q] indicates which 
values have been drawn fro m which source. Where avai l- 
able, we prefer the values of lLauberts fc Valent iin (1989), 
except in the case of ESO 323-G033, where, judging 
from inspection both of the Digitized Sky Survey plates 
and our acquisition i mage s, we consider the value of 
Ide Vaucouleurs et al.l (|1991l ) to be more appropriate. 

While i?25 is a reproducible measure of the radius of 
a galaxy disk, some caution must be exercised when 
comparing measurements in units of R25 between the 
present work and other studies of the outsk irts of galaxies 
(jFerguson et~allll998at fPohlen et alJl2002h . which gener- 
ally focus on face-on galaxies. In the absence of dust 
(an idealistic assumption, but more justified for the out- 
skirts of galaxies than for the gas-rich inner disk), sim- 
ple geometrical considerations lead us to expect the sur- 
face brightness along the line of sight to vary roughly 
as ~ l/cos(i), where i is the inclination angle. There- 
fore, a given isophotal radius will lie ~ 2.5 exponential 
scale lengths further out when the galaxy is seen edge- 
on than when it is seen face on. Because R25 typically 
corresponds to anywhere between 5 and 10 exponential 
scale lengths, R25 could change by 25 to 50% due to 
inclination. Therefore, quantitative comparisons of the 
distribution of outer-disk material in edge-on and face-on 
galaxies are complicated. 

The rotation curves for all the galaxies in our sam- 
ple are shown in Fig. OH] in units of the optical radius, 
i?25- We plot the measured line centroids for each indi- 
vidual row of the two-dimensional spectrum in which Ha 
is detected. We investigate every outlier and comment 
on the interesting ones below. Fluctuations of the rota- 
tion curves around the asymptotic velocity are typically 
larger than the uncertainties in the individual line cen- 
troid measurements. Statistical uncertainties in the cen- 
troid are below 10 km s _1 for regions with the strongest 
Ha flux, while the dispersion about the fitted curve is 
~ 20 km s _1 . It is possible that these small-scale fluctu- 
ations do represent actual kinemati c fluctuations across 
the s l it due to non-circular mo tions (jBeauvais fc Bothunl 
119991: lMcClure-Griffithil2006t ). However, it is also con- 
ceivable that these are due to surface brightness fluctua- 
tions across the slit. The slit image width of one arcsec is 
theoretically capable of introducing uncertainties in the 
line centroid of many tens of km s _1 , and it is there- 
fore not possible to analyze these small-scale and small- 
amplitude fluctuations any further, save to say that their 
velocity dispersion is an upper limit on the uncertainty 
introduced by surface brightness fluctuations across the 
slit. 

We have fitted functions of the form 

V(R) = V max x arctan{R/R )/(Tr/2) (1) 

to th e extracted rotation curves (|Stavelv-Smith et al.l 
using a simple \ 2 optimization with respect to 
the asymptotic rotation velocity V max , the scale length 
i?o, the central wavelength Ao, and the position of the 
galaxy center along the slit. Errors are a combination 
of the measured centroid uncertainties and the standard 
deviation of the data points around the fit in order to 
encompass the effect of surface brightness fluctuations 
across the slit and internal velocity fluctuations. This 
procedure implies that the overall quality of the fit can- 
not be judged from the \ 2 value (the reduced \ 2 wm 
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Fig. 2. — Velocity centroids of Ha measured relative to the sys- 
temic galaxy velocity on a row-by-row basis in the two-dimensional 
spectrum. The galaxies in this figure were observed with Gemini- 
S/GMOS. The solid line represents the best arctan fit. Open sym- 
bols show measurements where, to improve S/N, we have binned 
over larger regions of interest. 
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Fig. 3. — Velocity centroids of Ha measured relative to 
the systemic galaxy velocity on a row-by-row basis in the two- 
dimensional spectrum. The galaxies in this figure were observed 
with VLT/FORS1. The solid line represents the best arctan fit. 
Open symbols show measurements where, to improve S/N, we have 
binned over larger regions of interest. 



always be of order unity), but Figures [2] to [8] confirm 
that this functional form describes the shape of the in- 
dividual profiles well. Although this description is not 
a physical model, the fits provide two useful scaling pa- 
rameters for comparisons, V ma x and i?o- It is remarkable 
that an asymptotically rising (or, in the limit of small 
scale lengths, flat) rotation curve describes every object 
in our sample well; we do not see any significant evidence 
for declining rotation curves at large radii. This result 
may arise in part beca use our sample is dom inated by 
relatively faint objects (|Catinella et al .1 12006). but even 
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Fig. 4. — Velocity centroids of Ha measured relative to the sys- 
temic galaxy velocity on a row-by-row basis in the two-dimensional 
spectrum. The galaxies in this figure were observed with Magel- 
lan/B&C. The solid line represents the best arctan fit. Open sym- 
bols show measurements where, to improve S/N, we have binned 
over larger regions of interest. 
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Fig. 5. — Velocity centroids of Ha measured relative to 
the systemic galaxy velocity on a row-by-row basis in the two- 
dimensional spectrum. The galaxies in this figure were observed 
with MMT/Red Channel. The solid line represents the best arctan 
fit. Open symbols show measurements where, to improve S/N, we 
have binned over larger regions of interest. 



the brighter galaxies in our sample are consistent with 
an asymptotically rising or flat rotation curve out to the 
last measured data point. By visual inspection, only ESO 
445-G 081 offers indications of a turnover. The standard 
deviation of data points around the fit is on the order of 
~ 10 — 20 km s" 1 ; this contains contributions from the 
centroid measurement uncertainty as well, but is domi- 
nated by slit illumation effects and small-scale velocity 
fluctuations in the galaxy disk. 

Errors on t he best-fit p arameters are calculated using 
the method of lAvnil (|1976f ) by finding the upper and lower 
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Table 3. Rotation curve parameters 



Fig. 6. — Velocity centroids of Ha measured relative to the sys- 
temic galaxy velocity on a row-by-row basis in the two-dimensional 
spectrum for MCG -01-31-002, measured at Magellan with the 
B&C spectrograph. The solid line represents the best arctan fit. 
Open symbols show measurements where, to improve S/N, we have 
binned over larger regions of interest. 



limits to the error contours corresponding to a x 2 offset 
of 2.331 from the best-fit % 2 , appropriate for a fit with 
two degrees of freedom. We list the parameters resulting 
from the fit for each galaxy in Table [3] The asymp- 
totic velocities v max should not be taken as rotational 
velocities, especially in cases where our spectra have not 
probed the flat regime of the rotation curve (particularly 
ESO 380-G023). For purposes of recovering the rota- 
tional velocity, the fitting parameters should be used to 
calculate the velocity at a reasonably large radius, such 
as i?25, if a more uniform measure is desired. 

We do not observe significant correlations between the 
maximum extent of the Ha rotation curve and other 
galaxy parameters, such as i?-band luminosity or V C i rc . 

4.2. Notes on Individual Objects 

In this section, we briefly comment on specific obser- 
vations in individual galaxies. We only list objects here 
that deserve special attention due to peculiarities in their 
Ha rotation curve or other characteristics. A hard cut 
in signal/noise does not always perfectly separate plau- 
sible from spurious detections. Genuine outer-disk Ha 
emission can fall below a reasonable cut, while spurious 
surface brightness fluctuations can be found with very 
large signal/noise (e.g., as a result of background objects, 
CRs, or other image defects). It is therefore mandatory 
to inspect every spectrum visually to verify whether a 
weak detection (regardless of its kinematic properties) 
is plausible or not. While, in some of our spectra, we 
set a relatively low S/N threshold of 3 to flag possible 
detections, we adopt relatively conservative criteria for 
evaluating possible detections visually. In general, we 
consider a detection real if its line profile width is consis- 
tent with the slit image width, regardless of the position 
of the velocity centroid. 

4.2.1. ESO 323-G033 
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NOTE. — Parameters of the individual galaxy rotation curve, as 
gathered from the arctan fit. Columns show the asymptotic rotational 
velocity, the arctan scale length Rq in units of the isophotal radius 
-R25, and the central wavelength Ao- The sign of Rq only depends on 
the orientation of the slit relative to the sense of the rotation of the 
galaxy. 



This galaxy is one of the faintest (Mb ~ —17.7) and 
latest-type galaxies in our sample. It is also our deepest 
exposure. The rotation curve includes outer-disk emis- 
sion from at least five separate sources. Four of these are 
close to the disk (within i? 2 5, but distinct from the con- 
tinuous, bright, inner disk), while the outermost one is 
at a radius of ~ 1.5 x i?25- These outer-disk emission re- 
gions show striking kinematic anomalies: The outermost 
region is not consistent with the extrapolated rotational 
velocity of the galaxy at this radius, but rather is consis- 
tent with the systemic velocity of the galaxy itself. The 
innermost emission region is only significant at a level of 
2<7, and the formal uncertainties in the error centroid are 
large. Taken by itself, this data point would have to be 
disregarded, but the similarity of this detection to the 
four other emission regions, as well as its physical prox- 
imity to three of them, make it highly unlikely that this 
is spurious. The velocity centroid of this region is also 
displaced from the rotational velocity. 

To verify these detections and discuss whether the de- 
viation from the mean rotation curve fit could be caused 
by slit illumination effects, we show in Fig. [7] a detailed 
study of the kinematics of these five outer-disk emission 
regions. The bottom panel shows a greyscale projection 
of the relevant parts of the spectrum, including the Ha 
line as well as the neighouring [Nil] 6548 and [Nil] 6584 
lines visible as horizontal structures. Solid horizontal 
lines indicate the limits of the slit image at the rotation 
velocity extrapolated from the arctan fit. The upper pan- 
els show the derived velocity relative to the systemic ve- 
locity of the galaxy. The vertical structure in the bottom 
left panel is continuum light from a foreground star. 

All five emission regions are identifiable in this fig- 
ure. The outermost emission region is evidently dis- 
placed from the extrapolated rotation velocity because 
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Fig. 7. — Detailed study of Outer Disk Emission Regions in ESO 
323-G033. The upper panels show the rotation curve in two differ- 
ent radial regions of interest (the best-fit curve is plotted as a solid 
line), while the bottom panels show the reduced spectra. The pre- 
dicted position of slit-filling emission at the interpolated rotation 
velocity is indicated with two parallel, near- horizontal lines. Five 
distinct emission regions (the innermost one being very weak) are 
visible. The innermost and outermost regions deviate significantly 
from the extrapolated rotation velocity at their respective radii. 



the flux is almost entirely outside the region of the CCD 
that co-rotating gas in the slit would illuminate. For the 
innermost emission region, most of the flux also appears 
to fall well outside the relevant region. 

The integrated luminosity of the outermost emission 
region at a distance of 36.2 Mpc is ~ 4 x 10 35 erg s" 1 , 
about an order of ma gnitude faint er than M42, which 
has 4x 10 36 erg s" 1 (|Wilcotdll99l . The spatial extent 
along the slit is less than 2 arcsec, so that corrections for 
emission beyond the slit are unlikely to raise the estimate 
of the total luminosity by more than a factor of 2. This 
is considerably fainter than HII regions that are typically 
detectable by narrowband imaging of other galaxies, ow- 
ing to the fact that optical spectroscopy can achieve a 
much higher S/N per pixel than imaging for detecting 
line emission. 

The presence of gas near the systemic velocity at such 
large radii is surprising. One possibility is that we are ob- 
serving extreme cases of projection effects, i.e., gas that 
is not at its maximum elongation, so that only a small 
fraction of the orbital velocity is projected onto the line of 
sight. Given the close proximity and very homogeneous 
appearance of these emission regions in projection, it is 
unlikely that they should be so widely separated along 
our line of sight; furthermore, attempting to explain the 
low line-of-sight velocity of the outermost emission region 
by projection of a circular rotational velocity would re- 
quire the orbital radius to be so large as to be physically 
implausible (~ lli?2s)- We therefore reject the hypoth- 
esis that these emission regions are moving on circular 
orbits in the disk plane. While it is in principle possible 
that at least some of these outer-disk emission regions 
are satellite galaxies, which would allow them to be both 
at large radial distances and moving on non-circular or- 



bits, it is unlikely for even one satellite galaxy, let alone 
two or even five, to fall exactly on a one-arcsecond slit. 

There are two remaining plausible scenarios: First, the 
emission regions may be moving on strongly inclined or- 
bits, i.e., have a significant velocity component perpen- 
dicular to the plane of the galaxy, which would allow 
them to move on nearly circular orbits. However, this 
raises the question why five of these objects are seen 
within 1 arcsec of the plane of the galaxy, and why three 
of them have normal, disk-like kinematics. One possi- 
bility is that outer disk gas is colliding with material on 
inclined orbits (a galactic fountain or a recent accretion 
event may be a source of such material). In this case, 
shocks may be an important source of ionization lead- 
ing to the Ha emission. In fact, as will be discussed 
in a subsequent paper, the outer-disk emission-regions 
exhibit anomalously high [Nil] /Ha line ratios, which is 
consistent with the hypothesis that ionization mecha- 
nisms other t han local star formation are r esponsible for 
the emission (jBland-Hawthorn et al"]ll997f ). Second, the 
emission regions may be on highly eccentric orbits in the 
disk plane. In this case, they must have acquired their 
unusual kinematic properties within the past several hun- 
dred million years, as elliptical orbits would lead these 
emitters into the dense inner disk, where they would be 
destroyed quickly in collisions with the rotating gaseous 
disk. We may therefore be seeing infalling material. 

4.2.2. MCG -01-31-002 

This object is an optically inconspicuous, isolated Scd 
galaxy of M#=-19.6 at 5741 km s _1 , observed with the 
B&C spectrograph at Magellan. Data for the inner and 
outer disk come from two different exposures, the one 
in the outer disk being much deeper than the one of the 
inner disk, with a small gap in coverage just short of i?25- 
We detect gas at extremely large radii in this object, 
~ 2.5 x R 25 . 

In Figure [8] we show the Ha emission corrected to the 
mean fitted rotation curve along with cuts through the 
emission line at different galactic radii. The fluxes in the 
upper panels are scaled differently for display purposes, 
but the line flux in the three rightmost panels is more 
than an order of magnitude smaller than in the inner 
disk. The line profiles are fairly well-centered on the 
extrapolated rotation velocity and clearly following the 
disk rotation, with the possibility of a slightly higher- 
velocity tail at intermediate radii and a drop towards 
the outer edge. Most importantly, the line width of the 
faint Ha line in the outer disk is dramatically increased, 
compared to the bright Ha component in the inner disk. 

Could an instrumental misconfiguration, such as a 
wrong focus, be responsible for this result? The slit 
width, estimated from the sky lines, is 2.5 ± 0.1 pixels, 
while the Ha line in the central disk is 2.67 pixels wide. 
However, beyond i?25, we measure the width of Ha to be 
8.78 pixels, corresponding to a velocity broadening of 8.4 
pixels, 6.6 A, or ~ 300 km s _1 . The narrowness of the sky 
lines and the inner disk Ha suggest that the instrument 
is working fine. Despite the apparent high velocity dis- 
persion, the material is, on average, co-rotating with the 
disk. There is an exceptionally bright Ha knot (with a 
surface brightness approximately three times larger than 
anywhere else along the slit) at ~ 0.75 x i?25- Both the 
anomalously high velocity dispersion and the Ha hotspot 
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Fig. 8.— Detailed study of the spectrum of MCG 01-31-002. 
The left panel shows the spectrum relative to the fitted rotation 
curve. Data for the inner and outer disk come from two different 
exposures, the one in the outer disk being much deeper than the 
one of the inner disk, with a small gap in coverage just short of 
i?25- The right panels show the line profile in various bins, whose 
boundaries are indicated by the lines connecting the upper to the 
lower panels. The most important feature is a kinematically hot 
component of Ha emission, detected beyond i?25- Also notable is 
a bright Ha emission region near 0.75 i?25- 



may be signatures of on outer disk that was recently per- 
turbed and heated by tidal interactions or accretion. It 
is not possible to determine whether this broadened Ha 
component exists in the inner disk or is only present in 
the outer disk, as it is too faint to be visible on the short- 
term exposure centered on the galaxy center. 

4.2.3. UCG 10453 

This object, an Sd galaxy with Mb = —18.25, is 
remarkable for its very extended Ha disk (out to ~ 
1.7 x i?25)- The outer disk emission is clumpy, suggesting 
outer-disk star formation as the ionization mechanism. 
Only one side of the galaxy was sampled with deep ex- 
posures. 

4.2.4. ESO 380- G023 

At Mb ~ —17, this Scd galaxy at cz = 2, 753 km s _1 is 
among the faintest in our sample. The rotation curve is 
not very extended, but this object is notable for having 
the most-slowly rising rotation curve in our sample. Even 
at the largest radii to which we have sampled it (~ 0.9 x 
-R25), the rotation curve shows no sign of flattening. 

4.2.5. ESO 201-G022 

This Sbc galaxy with M B w -19 at cz = 4,068 km 
s , which is among our deepest exposures, shows Ha 
emission out beyond ~ I.3O-R25. It has a very regular 
rotation curve out to this radius. This galaxy is one 
of four targetted by Ide Griis et ail (|2001[) in a study of 
truncated disks. For this object, they find a truncation 
in the stellar surface brightness at ~ 100 — 110 arcsec 
with uncertainties on the order of 5-10 arcsec, which cor- 
responds to ~ 1.4i?25- Our outermost Ha detection is 



therefore consistent with being at a similar radius as the 
stellar truncation. 

4.2.6. ESO 299-G018 

The 3.4-hour exposure with FORS1 reveals emission 
extending at least out to 1.33 x i? 2 5, consistent with the 
extrapolated rotation curve, in this Sc galaxy of Mg = 
— 18. This object is problematic, as the Ha emission 
coincides with very faint sky lines, prompting us to set a 
relatively high detection threshold of la. 

4.2.7. ESO 381-G045 

This Sb galaxy is the earliest-type and highest-redshift 
object in our sample; at Mb = —19.6, it is also fairly 
luminous. The continuous rotation curve of the inner 
disk is not particularly extended; it reaches to 0.7 R25- 
However, the deep spectra reveal Ha (at 6 and 5cr, re- 
spectively) from two point-like emission regions that are 
clearly detached from the inner disk, extending the rota- 
tion curve to ~ 1.4. The luminosities (based on the flux 
within the slit) are slightly higher than that calculated 
for the outermost emission region in ESO 323-G033; they 
are 10 36 and 0.6 x 10 36 erg s . However, in this con- 
text, the reader is reminded again of the very uncertain 
flux calibration for the objects observed at Magellan; an 
uncertainty of a factor of 2 is possible in these estimates. 

4.3. Outer-Disk Detections 

The outermost detections of Ha emission and their 
kinematics are of particular interest. If they follow the 
mean rotation curve, then they set a lower limit on the 
extent of the disk, mark the material with (probably) 
the highest specific angular momentum, and provide the 
most extreme constraints on the rotation curve at large 
radii. To place the tightest constraints on these quanti- 
ties possible with our data, we use the following strategy: 
First, we determine the outermost easily detectable point 
of the Ha disk, but discarding isolated detections that are 
not continuous for at least several arcseconds. Where in 
doubt about the exact edge of the emission, we typically 
define larger spatial bins just before and beyond the pre- 
liminary edge, and measure fluxes and velocity centroids. 

We list in Table [4] the outermost Ha detections on ei- 
ther side of the galaxy separately, and include the mea- 
sured velocity relative to the interpolated rotation ve- 
locity at the bin center. In most cases, the outermost 
Ha detections are around i?25, although there are indi- 
vidual galaxies with significantly larger Ha disks, such 
as ESO 201-G022, ESO 299-G018, MCG -01-31-002, and 
UGC 10453, as well as galaxies with detached outer-disk 
emission regions (ESO 381-G045, ESO 323-G033). The 
velocity centroid uncertainties given in the table include 
contributions from surface brightness fluctuations across 
the slit, determined via the standard deviation around 
the fit as discussed before. Note that, for galaxies from 
the Magellan and MMT runs, the coverage is asymmet- 
ric, and deep exposures typically exist for only one side 
of the galaxy; the detections that are limited by radial 
coverage rather than by the significance threshold are 
marked. 

How well do these deep Ha rotation curves perform 
in comparison to optical spectroscopy with more conven- 
tional exposure times, as well as to 21-cm studies, and 
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are there indications that the true outer disks of galaxies 
are even larger than those observed here in Ha emission? 
To answer these questions, Fig. [3] shows the distribution 
of the outermost detections in our sample in comparison 

to several o ther sample s: 1) High-quality Ha rotation 

curves from lVogt et alJ (|2004l ). using exposures of up to 
3000 s duration on the Hale 200 inch telescope (itself 
a more powerful configuration than usual for the study 
of optical rotation curves) . Although predominantly fo- 
cused on clusters, this study also comprises a number of 
field galaxies, and we have selected only galaxies with 
morphological types later than Sb for the comparison. 
Instead of the -B-band R25 radius, this work uses the I- 
band R23.5 radius as a reference for the optical diameter. 

2) A lit erature compi lation of resolved 2 l-cm observa- 
tions bv lMartinl (|1998l ). We choose a subsample of spiral 
galaxies of Hubble type Sb or later at cz> 1,000 km s _1 
and require that the spatial resolution (i.e., the width 
of the half-power beam) be less than 20% of the opti- 
cal galaxy diameter. Although this sample is by its na- 
ture unsystematic, it provides an overview of the radii to 
which typical 21-cm observations probe outer disks. We 
use the optical radii listed in the Martin catalog itself, 
which are taken from the LEDA database and correspond 
to the the fj, = 25 mag arcsec" 2 isophotal surface bright- 
ness in the B-band, the same as adopted for our work. 

3) A 21-cm survey by iGarcfa-Ruiz. Sancisi fc Kuijkenl 
(2002) of edge-on late-type galaxies, aimed at the study 
of warps, and using the Westerbork Radio Synthesis Tele- 
scope. This sample is most similar to ours with regard 
to selection criteria. 

The optical rotation curves of lVogt et al.l (|2004f) gener- 
ally reach to .75 - 1.0 R25- The observations recorded in 
Martin (1998) typically detect HI to radii between 1 and 
1.5 times the optical diameter, although a few individual 
cases of extended hydro gen disks (~ 2 times the optical 
diame ter) exist. Finallv. [Garcia-Ruiz. Sancisi &: Kuijkenl 
(2002) typically detect 21-cm emission out to 1.5 R25, 
although a few cases of more extended emission exist. In 
comparison, of the deep optical spectra presented here, 
many can be traced out approximately to R25, but there 
is a significant number of galaxies for which Ha detec- 
tions extend considerably further, as far as 2 R25 in one 
case (out of 17), a region that was previously only probed 
in 21-cm; more than half of our sample show detections 
beyond R25. Therefore, in such cases, our observations 
are probing a regime that has previously not been ac- 
cessed routinely with optical spectroscopy. 

Figure [9] also shows that 21-cm observations on aver- 
age still probe even larger radii of the outer disk. This 
may have two reasons: First, Ha emission may sim- 
ply be too faint and drop below the detection thresh- 
old before the physical truncation of the neutral hydro- 
gen disk. Most ionization mechanisms (ionization by lo- 
cal star formation or b y UV flux from the inner disk 
(|Bland-Hawthornl [T998)) lead us to expect a declining 
surface brightness profile; only ionization by the cosmic 
ionizing background would eventually cause a flat sur- 
face brightness distribution as long as the disk is opti- 
cally thick to UV radiation, but there is no indication 
that our observations have reached the depth required 
for a detection of background- ionized gas. The second 
possible cause is that many, and possibly most, galax- 
ies are now thought to have warps in the outer disks 





i i i j i i i i | i i i i _ 
Vogt et al.(2004). Ha ~ 


1 t I .)■ r I 1 
1 1 1 1 1 1 kTSs^O 1 




1 1 1 1 1 1 1 1 1 1 1 
this work 1 


Martin (1998), 21-cm 

mi ^rWf^^lVi j 7 ? i i i - 


r Garcia-Ruiz et al. (2002), 21-cm: 



0.5 1 1.5 2 2.5 



R/R 26 

Fig. 9. — Comparison of the maximum extent (in units of R25) of 
our Ha detections (second panel from top) to Ha rotation curves 
with shorter exposure times from Vogt et al. (20(53) (top p anel), a 
literat ure compilation of resolved 21-cm observations from [Martini 
(1998) (third panel from top), and a 21-cm survey of lat e-type 
edge-on galaxies by Garcia- Ruiz, Sancisi & Kuiikcn (2002) (bot- 
tom panel). Our observations probe a domain normally not ob- 
served with optical spectroscopy; there is a significant subset of 
galaxies with Ha detections at radii previously only observed with 
radio telescopes. 



( Sanchez-Saavedra et al.l 120031: iReshetnikov et al. 


2002; 


Bosma 1982; IGarcfa-Ruiz. Sancisi & Kuiikenl 


2002), 



which means that, at large distances, the disks may phys- 
ically bend away from the plane observed by the spec- 
trographic slit. In fact, Christlein & Bland-Hawthorn 
(in prep.) are currently conducting a campaign to ob- 
serve warped galaxy disks with optical spectroscopy, and 
first obervations confirm that, in some cases, Ha emis- 
sion can be traced to larger distances in the warp than in 
the plane of the galaxy itself. However, we point out that 
the sample in the present study was specifically selected 
to exclude galaxies with indications of (optical) warps on 
publically available imaging. 

We find no obvious correlation between the maximum 
extent of the Ha detections and other galaxy properties 
such as absolute magnitude, morphology, or optical indi- 
cations of warps. 

What are the kinematic properties of the outermost 
Ha detections? The only detection that lags in velocity 
significantly from the fit is the outermost one in ESO 323- 
G033, which has been discussed in §4.2.11 There are no 
indications that the rotation curves turn over or decline 
out to the limits of the observable Ha distribution. This 
implies that at least out to these radii there is no strong 
truncation in the halo mass profile for any of these galax- 
ies and that the gas in the outer disk co-rotates with the 
inner disk and does not yet exhibit signs of transitioning 
to random, halo-like kinematics. ESO 323-G033 is the 
only exception to this second conclusion in our sample of 
17. 

4.4. Mass Density Profiles and Rotation Curves 

The existence of massive, extended dark matter ha- 
los around galaxies has been demonstrated using several 
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Table 4. Outermost Ha detections. 
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ESO 299- G 018 
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3.49 ±0.45 


1.18 


- 1.21 


-37.16 ± 18.00 


3.79 ±0.51 
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10.00 ± 15.19 


2.00 ±0.32 
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-8.38 ± 24.09 
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ESO 445- G 085 
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0.75 
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35.09 ±27.56 
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ESO 381- G 045 
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-9.79 ± 29.09 
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1.09 
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(0.64 - 
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41.40 ±9.26 2 


MCG -01-31-002 
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-13.66 ±24.40 
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(0.68 - 


0.70) 1 


19.75 ± 21.25 


57.14 ±8.80 2 


UGC 09138 
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UGC 09780 


1.03 - 
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32.61 ± 17.09 


9.13 ± 1.08 2 


UGC 10453 


1.57 - 


1.69 


36.73 ± 11.17 


1.24±0.35 2 
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3.84 ± 15.26 
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11.33 ± 13.07 


3.38±0.63 2 
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Note. — The first set of columns shows the outermost Ha detections for each galaxy with radius, velocity relative to the interpolated rotation 
velocity, and Ha surface brightness in 10~ 18 erg s -1 cm~ 2 arcsec - 2 . The second set shows the corresponding measurements on the opposite side 
of the galaxy. 1 This measurement is affected by the lack of radial coverage due to the short length of the slit, rather than the surface brightness 
level of the Ha emission. 2 Flux calibration for runs at MMT and Magellan is uncertain. 



lines of eviden ce (such as satellite galaxy kinematics an d 
weak lensing; IZaritskv et all Il993t iMcKav et al.l I2002T ). 
However, these are statistical studies, and it is there- 
fore difficult to exclude the possibility that there exists a 
small minority of galaxies wi th distinctly different prop- 
erties (jHonma fc Sofudll997f ). The kinematics of gas at 
large radii remains the most straightforward method of 
constraining the halo mass profile of individual objects 
because gas, if dynamically stable, must lie on closed or- 
bits. 

We pose two questions. Is there any case among the ob- 
served rotation curves that is consistent with a strongly 
truncated dark matter halo? Conversely, can a truncated 
halo explain the few instances of anomalous kinemat- 
ics in outer disk emission regions ijBland-H awthor n et al.l 
fl99l ? 

We attempt to answer these questions by first fitting 
a mass model that consists of an isothermal halo and a 
disk, to the observed rotation curves. The model has four 
free parameters, the halo central density and core radius, 
and the disk central surface density and scale length (the 
latter two could in principle also be constrained by the 
surface photometry, but we suspect that due to extinc- 
tion and stellar population gradients they are unreliable 
as indicators of the stellar surface mass density in this 
sample of nearly edge-on galaxies. Once we have a best- 
fit model, we fix these four parameters and introduce 
two additional parameters, the truncation radii of the 
disk and halo. Normalizing the probability for the un- 
truncated fit to unity, we reject cutoff radii for which the 
X 2 probability of the fit drops below 5%. 

We apply this analysis to the three galaxies with the 
most extended rotation curves: ESO 299-G018, UGC 
10453, MCG -01-31-002. We also include ESO 323-G033 
to verify whether a truncated halo could be responsible 
for the anomalous kinematic properties of at least the 
outermost Ha emission region. 

We present the results in Figure [10] For clarity, we 
show only a fraction of the measured rotation curve data 
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Fig. 10. — Disk-Halo models fitted to four of our most extended 
rotation curves. Solid lines show the best fit for untruncated disk- 
isothermal halo models. The dashed line and dotted line show 
models for a completely truncated halo, where the dashed line is 
a marginally acceptable model and the dotted line a marginally 
rejected one. In the best cases, outer-disk rotation curves rule out 
halo truncations within R.25- Even a complete truncation of the 
halo cannot explain the anomalously low velocity of the outermost 
Ha emission region in ESO 323-G033. 



points. We reach three conclusions: 1) truncations of 
the disk surface mass density have subtle effects on the 
rotation curve and are therefore not well-constrained, 2) 
at least in the most extended rotation curves in our sam- 
ple, dark matter halos cannot be truncated in any of 
these three galaxies at less than R ~ 0.9i?25, and 3) even 
a dark matter halo truncated well within i?25 has diffi- 
culties reproducing the unusual kinematic properties of 
the outermost emission in ESO 323-G033. 
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4.5. The Composite Spectrum 
4.5.1. Construction of the Composite Spectrum 

All rotation curves in our sample are described well by 
the same analytical form; their differences can be reduced 
simply to differences in the characteristic scale length, 
Ro, and the asymptotic rotational velocity V max . The 
homogeneity of the sample provides the opportunity to 
combine the spectra into a "composite spectrum" . Such 
a composite spectrum might allow us to detect finer and 
fainter details and characterize the behavior of an "aver- 
age" galaxy. 

In this section, we describe the construction of such a 
composite spectrum. We use the two characteristic scales 
obtained from the arctan fits to transform the individ- 
ual galaxy spectra onto a common coordinate frame be- 
fore superposing them. There are different choices avail- 
able for such a common coordinate system, and each one 
serves a particular purpose best. 

We use the radial projected distance from the galaxy 
center in units of R25 as the spatial axis. For the ve- 
locity axis, we use three mappings. First, for tracing 
the rotation curve to large radii and detecting emission 
with anomalous kinematics, we use the difference from 
the arctan fit, measured in km s _1 . This transformation 
is defined in such a manner that emission at rotational 
velocities greater than the interpolated rotation curve 
are plotted as positive and those which are less than the 
interpolated curve are plotted as negative. Second, we 
use the ratio of the observed velocity and the arctan fit. 
In this transformation, emission at the systemic velocity 
will always stack at an ordinate of y = 0, and counter- 
rotating gas will stack at an abscissa of y = — 1. This 
mapping is therefore ideal for searching for countcrrotat- 
ing gas. Third, to characterize the velocity dispersion of 
the gas, we determine the region illuminated in the two- 
dimensional spectrum by gas at the exact interpolated 
rotation velocity that fills the slit. We then assign for 
any detected emission beyond that region a y value that 
corresponds to the difference, in km s _1 , between the in- 
ferred velocity and the closest slit edge. This results in 
a minium possible velocity dispersion because it removes 
the maximal effect from uneven slit illumination. 

The flux uncertainties in the individual source pixels 
arc calculated as outlined in Sec. 13. 21 by determining the 
standard deviation of flux measurements carried out in a 
control region that is free of line emission or absorption. 
For determining the uncertainties in the centroid, we ap- 
ply a bootstrapping algorithm that randomly rcsamplcs 
the set of individual spectra from which the composite 
is constructed; for each resampling, the centroid is recal- 
culated, and the standard deviation after 100 iterations 
adopted as the centroid uncertainty. This estimate in- 
cludes uncertainties that are due to surface brightness 
fluctuations across the slit, and is therefore appropriate 
for evaluating whether there are radial trends of the ve- 
locity centroid in the composite spectrum. 

For the construction of the composite spectrum, we use 
all available spectra from the Gemini and VLT observ- 
ing runs, and the deep, outer-disk spectra from the MMT 
and Magellan runs. The fact that the latter do not pro- 
vide continuous spatial coverage gives rise to the concern 
that they could bias the shape of the surface brightness 
profile of the composite spectrum (a single galaxy with 
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Fig. 11. — Composite Spectrum and its surface brightness distri- 
bution. Ha emission is clearly discernible out to 1.5 X R25 in the 
bottom panel. Velocities in the bottom panel are are offsets from 
the predicted rotation velocity. The top panel shows the Ha sur- 
face brightness profile. Flux scale is linear below 1, and logarithmic 
otherwise (see text for detailed explanation) 



extremely low or high surface brightness in the MMT 
or Magellan samples could affect the outer parts of the 
composite spectrum and thus change the overall shape) . 
However, we have tested this possibility by calculating a 
composite spectrum from just the seven highest-quality 
VLT and Gemini spectra, and see no systematic differ- 
ence between both in either the overall shape of the spa- 
tial surface brightness distribution or the rotation curve, 
with the only exception that the Ha emission is more 
extended beyond 1.3 R25 in the composite of 17 spectra. 

4.5.2. The Extent of the Composite Spectrum 

The top panel of Fig. [TT] shows the Ha surface bright- 
ness as a function of radius. To include all of the data 
and make the low surface brightness levels legible, we 
plot values differently depending on the value itself as 
described by the following prescription: 

= f c(log(//c) + 1) if / > c; 
y \f/c if < c 

We set c = 1, i.e., fluxes up to 1 erg s _1 cm -2 arcsec -2 
will be plotted linearly, while every full step above 1 in- 
dicates an increase by a factor of 10. 

Ha emission is visible to a radius of at least 1.5i?25 in 
the composite spectrum. Between 1.25 R25 and 1.5 R25, 
Ha emission is significant at over 10 cr, but it drops to 
3(7 between 1.5 and 1.75 R25, which we do not accept 
as a significant detection. At smaller radii, there is no 
doubt that the emission up to 1.5i?25 is associated with 
the galaxy disks. Furthermore, the spectrum very clearly 
shows a broad absorption feature (indicated by the lower 
background level) surrounding the Ha emission within 
~ 0.75i?25- This feature is likely indicative of Ha ab- 
sorption. 

Individual galaxies exhibit a range of radii of their Ha 
disks, and in particular, emission in the composite spec- 
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trum at radii around 1.5i?25 is likely to be contributed by 
only a small subset of galaxies in the sample. Nonethe- 
less, even if we remove the galaxies with the most ex- 
tended individual Ha detections (ESO 201-G022, UGC 
10453, MCG -01-31-002, ESO 323-G033, ESO 299-G018) 
emission is still clearly discernible in the composite spec- 
trum out to ~ 1.25 x i?25- In general, though, emission 
in the composite spectrum is not more extended than 
in the most extended individual contributing spectrum, 
suggesting that faint outer-disk emission features are not 
a ubiquitous phenomenon whose S/N would improve as 
exposures are stacked. We reserve the discussion of the 
Ha surface brightness profiles, and in particular, how 
representative the composite is for the individual galax- 
ies, for a subsequent paper, save for the caveat that the 
kinematic constraints that we gain from the outermost 
detections in the composite spectrum are only obtained 
from a subset of galaxies. 

4.6. The Composite Rotation Curve 

We now examine the transformed two-dimensional 
spectrum in Figure [TT] for deviations from the fitted arc- 
tan rotation curve. Flux and centroid measurements are 
carried out in the same way as for the individual spectra, 
and uncertainties are calculated as described in M4.5.H 
The result is shown in Figure [T2] Detections below 3a 
significance are not plotted. 

Measurement uncertainties of the line centroid are 
shown in Figure [12] with a shaded area, and data points 
measuring the line centroid in larger bins of width 0.25 
R25 are overplotted. Out to 1.25 R25, the data points 
line up extremely well with AV^ = line, indicating that 
there are no systematic deviations from the arctan shape. 
The last significant data point between 1.25 and 1.5 R25 
exhibits a slight drop from the rotational velocity, but 
the uncertainties are large enough for this measurement 
to remain consistent with AV^ = 0. 

The extent of this average rotation curve does not pre- 
clude the possibility that individual galaxies have more 
extended outer disks, however, the composite spectrum 
demonstrates, better than most of the individual spec- 
tra, that rotation curves of late-type galaxies follow an 
asymptotically rising profile, at least to 1.5i?25, and that 
there are no significant systematic deviations from this 
profile in cither the positive or the negative direction. 
Nonetheless, it is important to be aware of a possible se- 
lection effect: the composite rotation curve at 1.5 R25 is 
contributed by only a small subset of galaxies; if kine- 
matic anomalies in the outer disk affect the maximum 
extent of the Ha emission (for example, due to warps), 
then such anomalies may be rendered undetectable for 
us. 

4.7. Velocity Dispersion in the Outer Disk 

Disk galaxies consist of two kinematically very distinct 
regimes — a disk dominated by ordered, circular mo- 
tions, and a halo dominated by random motions. Are 
these realms strictly separated from each other, or is 
there a transition region where formation mechanisms 
leave their signatures in the form of an increased veloc- 
ity dispersion? We have already seen possible examples 
of this interface region in individual spectra — the hot 
outer disk in MCG -01-31-002, and the gas on anoma- 
lous orbits in ESO 323-G0033. However, the question 
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Fig. 12. — Deviations of the Composite Spectrum from the arc- 
tan rotation curve. The bottom panel shows the rotation curve 
normalized composite spectrum. In the middle panel, we plot the 
velocity centroid of the Ha line in the composite spectrum relative 
to the arctan fit (solid line, surrounded by thinner lines marking 
the lcr standard deviation of the rotation curve relative to the fit). 
Velocity centroids have also been measured in several wide bins 
and overplotted as filled circles. Finally, the top panel presents the 
line profiles at several radii (each profile corresponds to one of the 
wide-bin data points in the middle panel). 



remains whether these are general features in the galaxy 
population, or whether they are rare and possibly short- 
lived phenomena that are only observable in a small frac- 
tion of objects. In this section, we addres this question 
by studying the velocity dispersion of the outer-disk gas 
along our line of sight in the composite spectrum. 

Optical spectroscopy has advantages and disadvan- 
tages for answering this particular question. Because the 
width of the slit is narrower than the typical beam size 
of 21-cm observations (even interfcrometric observations 
usually do not reach arcsecond-scale resolution), the tails 
of the observed line profile are representative of the veloc- 
ity dispersion rather than of a gradient in the line-of-sight 
velocity across the field (unless the galaxy is seen in an 
almost perfect edge-on configuration). Its drawback is, 
once again, the possibility that surface brightness fluc- 
tuations across the slit might be mistaken for velocity 
shifts. In a composite spectrum, these shifts might be 
mistaken for an increased velocity dispersion. To avoid 
this confusion, we use for this analysis the third type 
of mapping the individual spectra to the velocity axis 
described above, i.e., we cut out the region of each spec- 
trum that could be covered by emission anywhere in the 
slit moving at the rotational velocity. Any emission be- 
yond this region must originate from material that is not 
moving at the rotational velocity. 

In Figure [131 we show the newly constructed composite 
spectrum after removing the emission within the pro- 
jected slit image at the rotation velocity in each individ- 
ual spectrum. The top panel shows the ratio of Ha flux 
beyond the projected slit image relative to that within 
the slit image. The error bars represent two extreme 
cases: the upper error bar shows the ratio, assuming that 
the flux in the wings was underestimated by lcr, and the 
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Fig. 13. — The line-of-sight velocity dispersion of the disk gas, 
measured by the amount of flux beyond the predicted region of the 
slit image for gas co-rotating with the galaxy disk. The bottom 
panel shows the composite spectrum after the region of the line 
profile covered by the slit image has been removed from all indi- 
vidual spectra. Remaining flux emanates only from the wings of 
the observed line profile, which must be due to kinematic devia- 
tions from the fitted rotation curve. The upper panel shows the 
fraction of this flux compared to the total line flux. Solid lines 
and the shaded area show the ratio and its uncertainties measured 
on a row-by-row basis, while data points with error bars show the 
ratio measured in larger bins of width 0.25 R25. There is no solid 
evidence for an increase in the velocity dispersion with radius, al- 
though the final bin beyond 1.25 R25 is weakly constrained. 



total flux overestimated by la. The lower error bar is 
calculated analogously. 

The fraction of flux in the wings beyond the slit image 
is consistenly in the range between 0.2 and 0.3, with no 
indication of an increase up to a distance of 1.25 R.25- 
Beyond 1.25 R25, constraints are weaker; the last data 
point shows a rise, but is still consistent with the range 
seen in the inner disk due to the large uncertainties. Re- 
moving MCG -01-31-002, the galaxy with known high 
velocity dispersion, from the sample does not alter these 
results significantly; the recovered values shift within the 
measured uncertainties. We have also tested whether the 
flux in the line wings may be overestimated at large radii 
by the fact that the flux is always measured around the 
position of the strongest Ha detection; this carries the 
risk that the flux is evaluated at the position of a spuri- 
ous detection or at a position that is not even included 
in the original measurement of the full line flux. We 
have done so by a) measuring the flux at a fixed posi- 
tion in the composite spectrum, thereby eliminating the 
risk of "hunting for significance" , and by b) increasing 
the width of the measurement window that we use for 
detecting the full line flux. The measurements show no 
significant deviations of the recovered wing flux fractions 
from the presented case. 

Finally, the flux in the wings may be overestimated at 
small radii if the Ha emission follows the arctan profile 
better at some radii than at others. Systematic devia- 
tions from the arctan shape would manifest themselves 
by spilling more light beyond the calculated position of 



the slit image. To test this, we have constructed a new 
version of the composite spectrum in which the individ- 
ual spectra are not aligned along the interpolated rota- 
tion curve, but rather, along the empirically-determined 
line centroids. Repeating the analysis on this spectrum, 
the wing flux fraction changes by a small amount, but 
typically by no more than 0.05, which does not alter the 
observed trend. Only the data point between 1.25 and 
1.5 R25 shifts downward by almost 0.2, but this change 
is not significant. Therefore, this test yields no reason to 
change our previous conclusion that there is no increas- 
ing trend of velocity dispersion with radius. 

4.8. A Search for Counter- Rotating Gas 

We now turn from examining the bulk behaviour of 
the gas to examining whether there is a minority of gas 
with distinctly different kinematic properties. If galax- 
ies grow by hierarchical build-up, discrepant velocities 
should arise naturally as the galaxy aquires material with 
distinctly different angular momenta. Whereas in the 
dense, inner disk, such material is quickly assimilated 
through collisions with the gaseous disk, the lower den- 
sities and larger dynamic timescales at the outer edge 
of the disk will preserve such kinematic signatures. Ma- 
terial in counterrotation, in particular, provides striking 
evidence of hierarchical growth. 

Again, we use the composite spectrum to maximize our 
sensitivity to very faint Ha emission, and to that end, 
we transform the individual spectra into a coordinate 
frame where such emission is readily identifiable. For 
this purpose, we use the second mapping prescription, 
as described in M4.5.H by defining the coordinate along 
the vertical axis as the ratio of the measured velocity to 
the interpolated rotational velocity. In this coordinate 
system, material at the systemic velocity of the galaxy 
will always appear at a coordinate of 0, while counter- 
rotating material will be found at —1. To quantify the 
Ha emission on potentially counterrotating orbits, we in- 
tegrate the flux between —1.5 and —0.5 of the rotational 
velocity at each radius, and divide it by the flux between 
0.5 and 1.5 of the rotational velocity. 

In Figure [14] we plot the resulting ratio for the com- 
posite spectrum. At small radii, within ~ 0.4i?25, where 
the rotational velocities of the individual galaxies are still 
low, the tails of the co-rotating Ha emission line spill into 
the region where we search for counter-rotating gas, and 
therefore, these high ratios are not indicative of counter- 
rotation. At larger radii (beyond 0.5i?2s), the fraction 
of Ha emission with negative angular momentum is con- 
sistent with zero and suggests that no more than a few 
percent of the Ha emission comes from gas on such or- 
bits. Again, similar caveats apply as previously: at radii 
beyond 1.25 R25, few galaxies contribute to the com- 
posite, and if galaxies with counterrotating gas somehow 
have smaller Ha-bright disks, the true abundance of such 
kinematic features may be underestimated. 

5. CONCLUSIONS 

Based on the expectation that any accumulation of hy- 
drogen will emit Ha due to the omnipresence of ioniz- 
ing flux (due to local star formation, star formation in 
nearby galaxies, AGN, or the metagalactic UV flux), we 
obtained deep optical long-slit spectroscopic observations 
to detect Ha emission and to study the properties of the 
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Fig. 14. — Constraints on the mean amount of counter-rotating 
Ha-emitting gas in the composite spectrum. The bottom panel 
shows the composite spectrum, transformed into a coordinate sys- 
tem where gas at the circular velocity is at the ordinate +1, while 
counter-rotating gas is expected at -1. The upper panel shows the 
fraction of gas with counter-rotating kinematics, measured in the 
two strips around y = +1 and y = — 1 indicated in the bottom 
panel with horizontal lines. Within ~ 0.4 X i?25i the wings of the 
H« line profile of co-rotating gas are indistinguishable from gas on 
counter-rotating orbits, so that this figure should only be consid- 
ered at larger radii. Beyond 0.5 X R25, the fraction of Ha emission 
with negative specific angular momentum is small and consistent 
with zero. The solid line and shaded region show the measurements 
on a row-by-row basis, and data points with error bars are binned 
with a width of 0.25 R25. 



outer disks of galaxies. The purpose of this campaign 
was to characterize the kinematic behaviour of the outer 
disk gas and star forming regions, and in particular to 
search for kinematic anomalies in galaxies beyond the 
stellar disk that could be signatures of disturbances due 
to tidal interactions or mass accretion and might be the 
radial demarcation between galaxy disk and halo. We 
detect Hen emission down to fluxes of ~ several 10~ 19 
erg erg s _1 cm -2 arcscc~ 2 , and present in this paper 
our findings regarding the extent and kinematics of Ha- 
emitting gas in the outer disk. Further papers will use 
the same spectra to study the Ha and stellar surface 
brightness profiles as well as the metallicity of outer-disk 
gas. 

From our observations and analysis, we conclude for 
late-type, sub-L* galaxies as a class that: 



• Deep optical spectroscopy is a powerful comple- 
ment to 21-cm radio observations for studies of the 
outer gaseous disks of galaxies. Deep optical spec- 
troscopy offers an good combination of sensitivity, 
spatial resolution, and kinematic resolution, and, 
in the best cases, probes the baryonic disk to sim- 
ilar radii as typical 21-cm maps. In our composite 
spectrum, we trace Ha emission to 1.5 R25 (mea- 
sured edge-on) , and in individual galaxies, we reach 
~ 2i?25 in some cases, comparable to the largest 
radii traced in 21-cm studies. 

• There is a great variety in the distribution and ex- 
tent of Ha emission among outer disks. In at least 



one case, spatially continuous Ha emission extends 
to beyond 2i?25, while in others, outer-disk Ha 
emission is restricted to individual regions. Even 
in some of our deepest spectra, Ha is completely 
absent beyond R25. With our sample of 17 late- 
type objects, we are unable to relate this pattern to 
other galaxy properties, and it is unclear whether 
the lack of Ha emission at large radii in some galax- 
ies is due to a lack of ionizing flux or to an actual 
truncation or warp in the gaseous disk itself. 

• The kinematics of the outer disks are generally 
disk-like, with flat rotation curves and small ve- 
locity dispersion, out to the outermost edge of the 
detected Ha, typically up to 1.5i?25- The rotation 
curves in this sample rise asymptotically without 
major breaks or systematic deviations. 

• There is no evidence for an increase in the veloc- 
ity dispersion towards larger radii, as measured by 
the fraction of the Ha flux in the high- velocity tails 
of the Ha line. Since our survey is only sensitive 
to gas in the plane of the galaxy, deviations from 
the extrapolated rotation curve that may be asso- 
ciated with outer-disk warps may not be detected 
in this search; a follow-up program to study the 
kinematics of gas in warps is currently ongoing. 

• There is little, if any, counter-rotating gas at large 
radii. We place an upper limit of a few percent 
on such gas (under the prior that only gas that 
is in the plane of the galaxy and irradiated by a 
sufficiently strong ionizing flux is detectable to us). 
This result could mean that in these predominantly 
late- type galaxies tidal interactions or mergers that 
could give rise to counter-rotation are rare or that, 
even at these large radii, gas densities are suffi- 
ciently large to quickly eradicate gas moving on 
such anomalous orbits. 

In constrast to the general trends, we find several ob- 
jects that are notable for anomalous kinematics in the 
outer disk: 

• ESO 323-G033, a very late-type, faint object, ex- 
hibits several small, but homogeneous outer-disk 
Ha emission regions with luminosities as low as 
~ 10 36 erg s _1 . Two of these have velocities close to 
the systemic velocity of the galaxy, rather than the 
expected circular velocity at the respective radii. 
The similar appearance of all these emission re- 
gions, as well as their close proximity in projec- 
tion, make it unlikely that projection effects alone 
are responsible for the low line-of-sight velocities. 
It is more likely that they are either moving on ec- 
centric orbits in the plane, which implies that they 
cannot be old, or on orbits perpendicular to the 
plane, in which case it is likely that they are being 
highlighted as they impact gas in the outer disk. 

• MCG -01-31-002 exhibits an outer disk with an ex- 
ceptionally large velocity dispersion, on the order 
of the circular velocity, but which is nevertheless 
co-rotating. 
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These unusual kinematic featurs may join other phe- 
nomena, su ch as stellar stre ams observed in our own 
Milky W ay (iHelmi et al~ 20061) and around other, nearby 
galaxies ([Ferguson et al.T2002f ). as observational indica- 
tors of hierarchical growth. 

While many of the galaxies in our sample show regu- 
lar, disk-like kinematics as far as we can trace the Ha 
emission, the fact that at least 2 out of 17 objects show 
anomalous kinematics in the outer disk is particularly 
noteworthy, given that our targets were selected to be 
isolated late-type galaxies with no obvious signs of per- 
turbations or interactions. The true fraction of objects 
with such anomalies may lie even higher, as even our 
deepest spectra do not necessarily provide a complete 
census of the hydrogen in the outer disk. 

Based on observations obtained at the Gemini Ob- 
servatory (program ID GS-2005A-C-4), which is oper- 
ated by the Association of Universities for Research in 
Astronomy, Inc., under a cooperative agreement with 
the NSF on behalf of the Gemini partnership: the Na- 
tional Science Foundation (United States), the Particle 
Physics and Astronomy Research Council (United King- 
dom), the National Research Council (Canada), CON- 
ICYT (Chile), the Australian Research Council (Aus- 
tralia), CNPq (Brazil) and CONICET (Argentina). 
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